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Abstract 

We provide a physical basis for the local gravitational superenergy tensor. Further- 
more, our gravitoelectromagnetic deduction of the Bel-Debever- Robinson superenergy ten- 
sor permits the identification of the gravitational stress-energy tensor. This local gravi- 
tational analog of the Maxwell stress-energy tensor is illustrated for a plane gravitational 
wave. 
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The concept of energy and the law of conservation of energy play crucial roles in all 
physical theories. The notion of particle energy is already present in Newtonian mechan- 
ics. In more general physical theories, the notion of energy becomes more sophisticated, 
while still retaining its fundamental significance. In general relativity, the energy associ- 
ated with matter is represented by a stress-energy tensor T%f that allows us to define, in 
a consistent way, the local energy density of matter as measured by a given observer. The 
condition T^f.a = represents the law of local conservation of energy of matter, but it 
does not guarantee, in general, that the total energy is conserved. This is because T^f 
contains only the contribution of matter as well as all nongravitational fields, while the 
gravitational energy, which is also expected to contribute to the total energy, has not been 
taken into account. However, in general relativity there is no physically meaningful notion 
of energy of the gravitational field. In fact, one would expect that in general relativity the 
gravitational energy should be given by an expression quadratic in the first derivatives of 
the metric, such as the Landau-Lifshitz pseudotensor. However, it follows from Einstein's 
principle of equivalence that it is not possible to construct locally a covariant expression 
for energy only in terms of the spacetime metric and its first derivatives. In this paper, 
we present a physical derivation of an expression that is, in a certain average sense, the 
local stress-energy tensor of the gravitational field, an expression which is analogous to the 
stress-energy tensor of the electromagnetic field. In fact, the validity of this novel gravita- 
tional stress-energy tensor rests upon the approximate correspondence - in a quasi-inertial 
neighborhood surrounding the worldline of a geodesic observer - between gravitation and 
electromagnetism. Our gravitoelectromagnetic stress-energy tensor turns out to be directly 
proportional to the Bel-Debever- Robinson tensor, which is thereby furnished with a proper 
physical interpretation. 

The electromagnetic field is endowed with a symmetric traceless stress-energy tensor 
a/3 that is given in terms of the Faraday tensor F a/3 by 

where *F a p is the dual of F a p. In 1958, Bel suggested the notion of a gravitational 
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superenergy tensor by constructing an analogous fourth-rank tensor from the Riemann 
tensor R^ vpa , namely, 




(2) 



which is unambiguously defined only for spacetimes with = Ag pu [1]. In this paper, 
we assume that the cosmological constant vanishes ( A = 0) , and the signature of the Ricci 
flat spacetime is +2; moreover, units are chosen such that Newton's constant of gravitation 
and the speed of light in vacuum are set equal to unity. 

The gravitational superenergy tensor (2) has been used to define the local energy 
density and Poynting vector for source- free gravitational fields [1]. Contributions to this 
approach have also been made by Debever [2] and Robinson [3]; a detailed treatment of 
this topic as well as further references to the original literature is contained in Zakharov's 
monograph [4]. 

The Bel-Debever-Robinson (BDR) superenergy tensor T af3ry5 is totally symmetric and 
traceless. Moreover, it satisfies a conservation law T"^ 7 ^ = in analogy with 0°^ = for 
source-free electromagnetic fields. The BDR tensor has been particularly useful in numeri- 
cal relativity and has therefore been the subject of studies by a number of investigators [5] . 
However, a basic derivation as well as physical interpretation of the superenergy tensor has 
been lacking thus far. It is the purpose of the present work to ameliorate this situation. 

To provide a physical basis for the analogy with electrodynamics, essential use will 
be made here of the notion of a Fermi system. Consider an observer freely falling in a 
gravitational field; the observer carries a tetrad frame where A^ = dx^/dr is the 

observer's velocity and r is the proper time along its path. The Fermi system - which is 
the simplest generalization of a local inertial frame along the path of the observer - is a 
geodesic coordinate system based on a nonrotating frame along the observer's worldline 
[6]. Let X a = (t,X) be the Fermi coordinates along the observer's worldline; the spacetime 
metric in these coordinates is given by F g a f3 = Vap + F ^a/3 5 where rj^p is the Minkowski 
metric and F /i a ^(r,X) may be expressed as a power series in X away from the observer 
with coefficients that depend on the Riemannian curvature of spacetime. Only the lowest 
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order terms away from the path will be taken into account throughout this paper. It is 
useful to define the gravitoelectric and gravitomagnetic potentials $ g and A g , respectively, 
via F hoo = 2$ g and F hoi = —2(A g )i in analogy with electrodynamics. 
It follows from the construction of the Fermi system that 

d> 3 = -\ F R 0i0j (r)X i X\ (3) 

(A g )i = ± F R 0j i k (r)X^X k , (4) 

where F R ce ^s(r) = R(a)(0)(>r)(6) = R ^P^( a )^((3)^)^(5) is the curvature as measured by 
the fiducial observer. A discussion of the general properties of this curvature is beyond the 
scope of this paper [7,8]. We now define the gravitoelectric {£) and gravitomagnetic (£>) 
fields in complete analogy with electrodynamics; hence, 

£i(r,X)= F R 0i0j (T)X\ (5) 

Bi(r,X) = ~ e ijk F R jk0 i(r)X l , (6) 

which agree with the identification of "electric" and "magnetic" components of the Rie- 
mann tensor. The curvature tensor generally consists of "electric" , "magnetic," and "spa- 
tial" components; in a Ricci flat region, however, the "spatial" components are given by the 
"electric" components. Therefore, Eqs. (5) and (6) contain the full information regarding 
the gravitational field along the path. Moreover, it is possible to express the gravitational 
field equations in the Fermi system in a form that is analogous to the Maxwell equations 
using Eqs. (5) and (6). This partial agreement between the field equations of gravitation 
and electrodynamics is well known and will not be further elaborated here [9]. 

The Lorentz force law is needed in addition to Maxwell's equations to produce a 
complete picture of classical electrodynamics. In the Fermi system, this equation to first 
order in velocity is 
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where m is the inertial mass of a test particle, and q E and q B are the gravitoelectric 
and gravitomagnetic charges of the particle, respectively. It follows from the gravitational 
Larmor theorem [10] that q E = —m and q B = —2m; the negative signs of the gravito- 
electromagnetic charges reflect the attractive nature of gravitation, while q B /q E = 2 since 
gravitation is a spin-2 field. Equation (7) turns out to be the generalized Jacobi equa- 
tion [11], which is the equation of motion of a free test particle relative to the fiducial 
observer, valid to first order in X and (iX/c/r. This agreement between the deviation 
equation and the Lorentz force law is a remarkable result, since it makes it possible, in 
principle, to introduce the stresss-energy tensor for gravitoelectromagnetism in complete 
analogy with the standard deduction of the electromagnetic stress-energy tensor. However, 
the electromagnetic stress-energy tensor is defined globally while the gravitoelectromag- 
netic stress-energy tensor would have to be defined in a narrow tube along the timelike 
path of the observer in spacetime. The invariance under translations in time and space 
is ultimately responsible for the existence of conservation laws of energy and momentum, 
respectively, as well as the stress-energy tensor in Minkowski spacetime. The existence of 
the gravitational pseudotensor can be similarly justified in an asymptotically Minkowskian 
spacetime manifold. On the other hand, the gravitoelectromagnetic stress-energy tensor 
would owe its existence to an inertial Fermi region (called the "tidal frame" in [12]) in the 
form of a thin Minkowskian cylinder in spacetime along the timelike path of the observer. 
It is possible to define a gravitational Faraday tensor in the Fermi system by 

Fafi = — F RapoiX , (8) 

which contains Eqs. (5) and (6). Let us define the stress-energy tensor associated with 
T a f3 in the standard manner by 

4tt 4 
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which is symmetric and traceless by construction. The substitution of Eq. (8) in Eq. (9) 
reveals that T af3 is zero along the path of the fiducial observer (X= 0) ; this is ultimately a 
consequence of Einstein's principle of equivalence. On the other hand, T a ^ is nonzero in a 
neighborhood of this path as the gravitoelectromagnetic field initially varies linearly with 
distance away from the fiducial observer. Imagine a second observer with a worldline only 
infinitesimally distant from the fiducial observer; if a second Fermi system is established 
along this path and our construction of the gravitoelectromagnetic stress-energy tensor is 
repeated, then this tensor would vanish identically along the second path but would give 
a finite result along the original reference path. It follows from these considerations that 
the proper stress-energy tensor along the path of an observer must be defined through an 
averaging procedure. At any given proper time r, let us average the tensor given by Eq. 
(9) over a limiting tube along the worldline. More precisely, consider a sphere of proper 
radius eL, 0<e<l, centered around X= at a given r such that 

T (a)(/3) (r) = Um £ -2 < T a/3 > (1Q) 

where the angular brackets denote the operation of averaging over the sphere. Here L is 
a constant invariant length scale that is otherwise unspecified; in practice, it could be a 
characteristic length scale for the problem under consideration (e.g., in the exterior field of 
a black hole, L could be the mass of the black hole), or, in the absence of a natural scale 
in spacetime as in a pure gravitational radiation field, one could set L equal to the Planck 
length. Using the fact that < X l X^ >= (e 2 L 2 /3)5ij, we can express Eq. (10) in arbitrary 
coordinates as 

L 2 

= Y^ T M^ A (0) A (0)> (11) 

where T^ vpcj is symmetric and traceless in the first pair of indices and symmetric in the 
second pair of indices by construction, and is given by 
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This tensor is identical with that first defined by Bel [1] in 1958 for = Ag^ and coin- 
cides with Eq. (2). It is important to remark, however, that in our gravitoelectromagnetic 
deduction, i.e. Eqs. (3) - (12), no restriction has been placed on the Ricci tensor; therefore, 
Eqs. (11) and (12) define the gravitoelectromagnetic part of the local stress-energy tensor 
of a general gravitational field. Furthermore, the linear treatment of the fields in Eqs. 
(5) and (6) is sufficient to obtain the general results in Eqs. (11) and (12) since terms of 
higher order would simply drop out of Eq. (10) in the limit e — > 0. 

To define the local gravitational stress-energy tensor at an event in spacetime, a 
geodesic observer with ideal gyroscope axes is needed at the event; then, T^ a )(p) = 
L 2 T( Q ,)( (S )(o)(o)/l27r for the observer, and the stress-energy tensor for any other observer at 
that event can be obtained from by an appropriate Lorentz transformation. More- 

over, the local physical quantities defined in this way (i.e. energy density, energy flux, 
momentum density, and the gravitational stresses) have now their proper dimensionality 
as a consequence of the introduction of the constant length L. The freedom in the choice 
of L implies, however, that these quantities are fixed up to a constant scale factor. 

To illustrate these results, let us consider a plane gravitational wave given by 

-ds 2 = -dt 2 + dx 2 + U 2 (e 2h dy 2 + e~ 2h dz 2 ), (13) 

which propagates along the x-axis and is linearly polarized. Here U and h are functions 
of it = t — x and are related by U" + h' 2 U = 0, where a prime indicates differentiation 
with respect to u. This gravitational field is of Petrov type N and is a member of a class 
of plane wave spacetimes that has been extensively studied [13]. 

In the coordinate system under consideration here, observers located at fixed spatial 
coordinates follow geodesies. Furthermore, the natural tetrad system A^ associated with 
these observers is nonrotating. The nonzero elements of are given by A^ = A^ = 
1, X 2 2 ^ = U~ l exp(— h), and A^ = U~ l exp(/i). The local gravitational stress-energy 
tensor according to these observers is given by 
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K{u) = h" + 2^-ti. (15) 



The spacetime is flat if K{u) = 0. All the nonzero components of ^(a)(/3)(7)(*) can be 



obtained from Eqs. (14) and (15). As expected, the only nonzero component of the 
gravitational stresses is the local pressure of the radiation along the x-axis. 

Let p g be the local energy density of the gravitational radiation field according to the 
fiducial static observers; then, p g = L 2 K 2 (u)/6tt. What is the energy density p g measured 
by observers boosted along the x-axis with speed /?? The boost can be expressed in 
terms of coordinates as u = Au, Av = v, y = y, and z = z, where v = t + x and 
A = (1 — (3) 1 ' 2 /{I + ff) 1 ! 2 is the Doppler factor. By transforming the tetrads at a given 
event, one can show explicitly that p g = A 4 p g . Alternatively, this result may be obtained 
by noting that under the boost the metric (13) and hence the magnitudes of U and h 
remain invariant; then, the form of Eq. (15) leads directly to the formula for p g . 

It follows that the local gravitational radiation energy density transforms in this case 
like the square of the energy density in standard field theory in Minkowski spacetime. This 
could possibly have observable consequences for the local inertia of gravitational radiation; 
however, such radiation can not be physically confined as a consequence of the universality 
of gravitational interaction. 

The application of these concepts to rotating gravitational waves [14] will be published 
elsewhere. 

Two of the authors (B. M. and J. C. M.) gratefully acknowledge the hospitality ex- 
tended to them at ICN-UNAM, Mexico City, where part of this work was done. We wish 
to thank R. Sussman for helpful discussions. 



8 



References 



[1] L. Bel, C. R. Acad. Sci. Paris 247 (1958) 1094; 248 (1959) 1297; Cah. Phys. 16 
(1962) 59. 

[2] R. Debever, Bull. Soc. Math. Belg. 10 (1958) 112. 

[3] See the paper of L. Bel, La Radiation Gravitationnelle, in: Colloques Internationaux 
du CNRS, Les Theories Relativistes de la Gravitation (CNRS, Paris, 1962) pp. 119- 
126. 

[4] V. D. Zakharov, Gravitational Waves in Einstein's Theory, translated by R. N. Sen 

from the 1972 Russian edition (Halsted, New York, 1973). 
[5] J. A. Wheeler, Phys. Rev. D 16 (1977) 3384; A. Peres, Phys. Rev. D 18 (1978) 608; 

J. Ibanez and E. Verdaguer, Phys. Rev. D 31 (1985) 251; J. Isenberg, M. Jackson and 

V. Moncrief, J. Math. Phys. 31 (1990) 517; N. Breton, A. Feinstein and J. Ibanez, 

Gen. Rel. Grav. 25 (1993) 267. 
[6] For recent interesting discussions of the Fermi system see K.-P. Marzlin, Phys. Rev. 

D 50 (1994) 888; Gen. Rel. Grav. 26 (1994) 619. 
[7] J. K. Beem, P. E. Ehrlich and K. L. Easley, Global Lorentzian Geometry, second 

edition (Marcel Dekker, New York, 1996) ch. 2. 
[8] B. Mashhoon, Phys. Lett. A 163 (1992) 7; B. Mashhoon and J. C. McClune, Mon. 

Not. Roy. Astron. Soc. 262 (1993) 881. 
[9] A. Matte, Canadian J. Math. 5 (1953) 1; B. Mashhoon, Found. Phys. 15 (1985) 497; 

R. T. Jantzen, P. Carini and D. Bini, Ann. Phys. (NY) 215 (1992) 1; W. B. Bonnor, 

Class. Quantum Grav. 12 (1995) 1483. 
[10] B. Mashhoon, Phys. Lett. A 173 (1993) 347. 

[11] B. Mashhoon, Astrophys. J. 197 (1975) 705, N.B. Eq.(A7) of the Appendix; D. E. 

Hodgkinson, Gen. Rel. Grav. 3 (1972) 351. 
[12] B. Mashhoon, Astrophys. J. 216 (1977) 591. 

[13] H. Bondi, F.A.E. Pirani and I. Robinson, Proc. Roy. Soc. London A 251 (1959) 519. 
[14] B. Mashhoon and H. Quevedo, Phys. Lett. A 151 (1990) 464; H. Quevedo and B. 
Mashhoon, in: Relativity and Gravitation: Classical and Quantum, Proc. SILARG 

9 



VII, eds. J. C. D'Olivo et al. (World Scientific, Singapore, 1991) pp. 287-293. 



10 



